Effective convergence of the 2PI-l/iV expansion for nonequilibrium quantum fields 
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The 1/N expansion of the two-particle irreducible effective action offers a powerful approach to 
study quantum field dynamics far from equilibrium. We investigate the effective convergence of the 
1/N expansion in the O(N) model by comparing results obtained numerically in 1 + 1 dimensions 
at leading, next-to- leading and next-to-next-to- leading order in 1/N as well as in the weak coupling 
limit. A comparison in classical statistical field theory, where exact numerical results are available, 
is made as well. We focus on early-time dynamics and quasi-particle properties far from equilibrium 
and observe rapid effective convergence already for moderate values of 1/N or the coupling. 

PACS numbers: 



Introduction - The need to understand the dynamics 
in heavy ion physics and the quark gluon plasma as well 
as highly nonlinear phenomena in the (post-)inflationary 
universe has lead to substantial development in the study 
of quantum field evolution far from equilibrium. One par- 
ticular approach, firmly based on functional methods in 
field theory, employs the two-particle irreducible (2PI) 
effective action The nonperturbative 2PI-1/JV ex- 
pansion 0, where N denotes the number of matter 
field components, has in particular proven fruitful. In 
recent years this approach has been applied to a variety 
of nonequilibrium problems in 3 + 1 dimensions, related 
to inflationary preheating [4L a], effective prethermaliza- 
tion fermion dynamics [7 1, transport coefficients [H 
and kinetic theory 0, [1(3] , slow-roll dynami cs fll| , topo- 
logical defects [12J. nonthermal fixed points [13j , expand- 
ing backgrounds [14} . nonrelativistic cold atoms |15j, etc. 
Besides these applications, theoretical aspects are under 
continuous investigation: renormalization in equilibrium 
is by now well understood 1|| 17, 18 1 and practical imple- 
mentations of renormalization out of equilibrium are be- 
ing developed [l9(. Moreover, progress in adapting these 



methods to gauge theories is steady [2CJ, |21|, [2 



Applications of 2PI effective action techniques are 
based on truncations, employing either a weak coupling 
or a large N expansion. So far, truncations stop at rel- 
atively low order: as far as we are aware all studies in 
field theory use next-to- leading order (NLO) truncations 
in the coupling or 1 /N [3fJ . The reason for this situation 
is clear: beyond NLO the complexity of the equations 
and the numerical effort required to solve them increases 
dramatically [23j. In this paper we present the first re- 
sults beyond NLO in field theory. 

Remarkably, the lowest order truncations beyond mean 
field theory include already many of the physical pro- 
cesses necessary to describe quantum field dynamics both 
far and close to equilibrium and are capable of capturing 
effective memory loss, universality of late time evolution, 



and thermalization. The natural question to ask is then 
how accurate a truncation at a given order describes the 
dynamics in the full theory. There are several ways this 
can be investigated. In the case of a systematic expan- 
sion, it should be possible to compare different orders 
of the expansion, shedding light on the effective conver- 
gence. When restricting to LO and NLO truncations 
only, the applicability of this approach is limited. The 
reason is that in LO mean field approximations scattering 
is absent and there is no notion of equilibration and ther- 
malization, as there is at NLO. It is therefore necessary to 
consider the next-to-next-to-leading order (N 2 LO) con- 
tribution as well. In the first part of this paper we study 
this problem in the O(N) model and compare the dy- 
namics obtained at LO, NLO and (part of) N 2 LO in the 
1/N expansion in quantum field theory. In cases where 
an exact solution is available, a direct comparison can 
be carried out, and this approach has been successfully 
applied using classical statistical field theory instead of 
quantum field theory 25|, 26}. In the classical limit, the 
nonperturbative solution can be constructed numerically 
by direct integration of the field equations of motion, 
sampling initial conditions from a given probability dis- 
tribution. We use this approach to further quantify the 
role of truncations in the second part of the paper. 

2PI dynamics - We consider a real TV-component scalar 
quantum field with a A(<^> a </) a ) 2 /(4!iV) interaction (a — 
1, . . . , N). In order to allow for an economical 1/N ex- 
pansion, an auxiliary field x is introduced to split the 
four-point interaction, which results in the action 



S[(f>,x] = - J dx° 



rfx 



% N 2 1 , - ' 
~2XX + -jXVa^a 



(1) 



The theory is formulated along the Schwinger-Keldysh 
contour C, as is appropriate for initial value problems. In 
the symmetric phase ((<f> a ) = 0), the 2PI effective action 



depends on the one-point function \ — (x) and the two- 
point functions 

G a b(x,y) = G(x,y)S ab = (T c <pa(x)4>b(y)) , (2) 
D(x,y) = (T cX (x)x(y)) ~ (X(x))(x(y))- (3) 
In terms of those, the action is parametrized as [H, [27| 

T[G,D, X ] = S[0,x] + ^TrlnG- 1 + ^TrG^ 1 (G - G ) 



--TrlnD- 1 
2 



-TrD^(D-Do) 



T 2 [G,D], (4) 



where G 1 = i(D + m 2 + x) and D 1 = 3N/(i\) denote 
the free inverse propagators. Extremizing the effective 
action gives equations for x, G and D. The latter take 
the standard form G -1 = G 1 - E and I?" 1 = Z^ 1 - II, 
with self energies E = 2i(5F 2 /(5G and II = 2iST 2 /SD. 

In order to solve an initial value problem, the prop- 
agators and self energies are decomposed in statistical 
(F) and spectral (p) components. For instance, the ba- 
sic two-point function is written as G(x,y) — F(x,y) — 
(i/2)sign c (a; — y°)p(x,y), and the causal equations for 
F and p take the form 



p x +M 2 (x)] F(x,y) 



dz E p (x, z)F(z,y) 



dzT, F (x,z)p{z,y), (5) 



p x + M 2 (x)]p(x,y) = - / dzE p 

JyO 



{x,z)p{z,y). 



Here we used the notation f x dz — ( x dz° [ dz. The 

Jy Jy J 

effective mass parameter is given by M 2 (x) = m 2 + 
X = m 2 + \(N + 2)/(6N)F(x,x). After decomposing 
D (x, y) = X/(3N) [iS c {x-y) + D F (x, y) - (i/2)sign c (x° - 
y°)D p (x, y)] , similar equations can be found for Dp and 
D P I- 

2PI-1/N expansion - In the 2PI-1/AT expansion the 
nontrivial contribution T 2 [G, D] is written as T 2 = 
pNLO _|_ pN lo _j_ _ _ Powercounting is straightforward: a 
closed loop of G propagators yields a factor N, whereas a 
D propagator contributes l/N. This yields one diagram 
at NLO (~ iV°) and two diagrams at N 2 LO (~ l/N), see 
Fig.ffl Cutting &G/D line gives the self energy E/II. At 
NLO the self energies have no internal vertices and are 
therefore easily evaluated in real space, where they are 
given by the product of two propagators. For example, 
the statistical component of the NLO self energy is given 
by the expression 



: N F LO (x,y) = -g 



F(x, y)D F {x, y) - -p(x, y)D p (x, y) 



(6) 

where g — A/(37V). At N 2 LO, the number of terms in- 
creases substantially, due to the possibility that any line 






(2) 



(3) 



(4) 



FIG. 1: 2PI-1/JV expansion: NLO (2 loops) and N 2 LO (3 
and 4 loops) contributions. The full/dashed line denotes the 
G/D propagator. 



can be of the F or p type. Moreover, self energies have 
two or four internal vertices, greatly increasing the com- 
plexity of the expressions. Let us first consider the 3-loop 
diagram at N 2 LO. The corresponding self energy has two 
internal vertices and five propagators. The complete ex- 
pression can be found in Ref. [23j and has O(30) terms. 
Here we give two terms to illustrate the structure, 

E^ LO (x,y) = -g 2 I dz I dw p(x, z)D F (x, w)p(y,w) 



D F (y, z)F(z, w) + -D p (y, z)p(z, w) 



+ 



(7) 



The two internal vertices result in two nested memory 
integrals. The second diagram at N 2 LO yields self ener- 
gies with four nested memory integrals. In contrast to 
NLO, at N 2 LO the evaluation of self energies completely 
dominates the numerical effort. 

Quantum dynamics far from equilibrium - We con- 
sider a spatially homogeneous system and solve the dy- 
namical equations numerically in 1 + 1 dimensions, by 
discretizing the system on a lattice with spatial lat- 
tice spacing a and temporal lattice spacing a t , using 
at /a — 0.2. Initial conditions are determined by a 
Gaussian density matrix, resulting in initial correla- 
tion functions F(0, 0;k) = F(at,a,t;k) = [no(wk) + 
l/2]/wk, F(a t ,0;k) = F(0, 0; k)(l - a 2 u 2 /2). Here 
no is the initial (nonequilibrium) particle number and 

1 /2 

cjk = (k 2 + Mq) with Mq the initial mass, determined 
sclfconsistently from the mean field mass gap equation. 
Initial conditions for the spectral function arc determined 
by the equal time commutation relations. 

We consider LO mean field dynamics, obtained by 
putting the nonlocal memory integrals on the right-hand 
side of Eq. ([5|) to zero. The only nonequilibrium aspect 
is in the time dependent mass M 2 (x). The NLO ap- 
proximation is solved without further approximation. At 
N 2 LO the numerical effort differs substantially between 
the 3-loop and the 4- loop diagram in Fig. [TJ due to 2 and 
4 internal nested memory integrals respectively. There- 
fore we restrict ourselves to the first diagram at N 2 LO 
and refer to this as N 2 LO'. In principle, a subtle can- 
cellation between the 3- and 4-loop diagram could occur. 
However, since the second N 2 LO diagram is (naively) 
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FIG. 2: Unequal time correlation function F(t, 0; k = 0) for 
N = 2,4, 10, 20 in the quantum theory. The dashed/full lines 
show results for the 2PI-1/W expansion at NLO/N 2 LO'. 
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FIG. 3: Masses and widths extracted from the unequal-time 
correlation functions in Fig. |5|as a function of 1/N. In the 
mean field approximation, the width is zero for all values of N. 



suppressed by two powers of the coupling constant, this 
is found not to be the case. We come back to this below. 
A similar comparison between 2- and 3-loop truncations, 
appearing at the same order in a coupling expansion, was 
made in Ref. 



24] 



The unequal-time two-point function at zero momen- 
tum is shown in Fig. [5] The coupling is A/m 2 = 30, 
where m is the renormalized mass in vacuum. The lattice 
spacing is am = 0.5 and volume Lm = 32. The mem- 
ory kernel is preserved completely. At LO (not shown) 
there is no damping in the unequal-time correlation func- 
tion. Beyond LO, we observe an underdamped oscilla- 
tion, the signal of a well-defined quasiparticle. Increasing 
N shows convergence of the two truncations considered. 
We can extract the quasiparticle properties far from equi- 
librium by fitting these curves to an Ansatz of the form 
Ae~ 7t cos Mt, where M is the quasiparticle mass and 7 
its width (divided by 2). The resulting values are shown 
in Fig. [3] as a function of 1/N. We observe effective con- 
vergence when the value of N is increased, as expected 
for a controlled expansion. For N > 10, the NLO and 
N 2 LO' results are practically indistinguishable. In 1 + 1 
dimensions the perturbative onshell width from 2^2 
scattering processes is not defined due to the constraints 
from energy-momentum conservation. However, includ- 
ing higher order effects via the use of selfconsistently de- 
termined propagators leads to a width that vanishes as 
1/N, as expected from naive powercounting. The quasi- 
particle mass converges to the value determined by the 
nonthermal fixed point of the mean field equations in the 
large N limit, which can be computed analytically [25j |. 
For the parameters used here we find M/m = 1.48 when 
N -> 00. 

Classical statistical field theory - In order to further 
investigate the applicability of the 1 /N expansion and in 
particular the role of the 4-loop diagram at N 2 LO, we 



turn to classical statistical field theory, where the non- 
perturbative evolution can be computed by numerically 
solving the classical field equations, sampling initial con- 
ditions from the Gaussian ensemble. In order to carry out 
the comparison, we also take the classical limit in the set 
of 2PI equations. In this limit statistical (F) two-point 
functions dominate with respect to spectral (p) functions. 
There are several ways this can be motivated, using for 
instance classical statistical diagrams or arguing for clas- 
sicality when occupation numbers are large 



26, 28, 2 



The result is that terms are dropped that are subleading 
when F is taken to be much larger than p. For example, 
both in Eq. © at NLO and in Eq. © at N 2 LO, the last 
terms are dropped with respect to the first ones. We have 
carried out this procedure for all terms appearing in self 
energies at N 2 LO' in Ref. 

Since the 4-loop diagram at N 2 LO is suppressed by 
the coupling constant, its effect is expected to diminish 
at weaker coupling. In order to verify this, we use a rela- 
tively small value of N = 4, where NLO and N 2 LO' differ 
at large coupling, both in the quantum and the classical 
theory. A comparison between NLO, N 2 LO', and the ex- 
act numerical result (MC for Monte Carlo) is shown in 
Fig. |U In the numerical integration of the classical equa- 
tions of motion, we have used a sample of 2.5 x 10 5 initial 
conditions. In the top left corner, the coupling constant 
is A/m 2 = 30. Damping at N 2 LO' is less than at NLO, 
similar to the situation in the quantum theory. Decreas- 
ing the coupling constant, we observe that the agreement 
between the two 2PI truncations improves. Moreover, for 
weaker coupling, we note that the MC results lie between 
the curves at NLO and N 2 LO', signalling a well converg- 
ing expansion. For even smaller coupling (not shown), 
we found that all curves lie on top of each other, even for 
N = 4. Since the N 2 LO approximation is numerically 
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FIG. 4: Unequal time correlation function F(t, 0; k = 0) for 
N = 4 in the classical theory. Included are the NLO and 
N 2 LO' results in the classical limit of the 2PI-l/jV expansion, 
and the exact result obtained by numerical integration of the 
classical equations of motion (MC). The coupling constant is 
A/m 2 = 30/n, with n = 1, 2, 4, 8. 



substantially more intensive than the NLO one, we have 
investigated whether the memory kernel at N 2 LO could 
be truncated earlier than at NLO, saving considerable 
computer resources. However, we found this not to be 
the case: in the results shown here it was necessary to 
preserve memory kernels over the entire history, also for 
the N 2 LO contribution. At this moment this prevents us 
from going to larger times and studying the approach to 
equilibrium beyond NLO. 

Summary - We have investigated convergence proper- 
ties of the 2PI-1/7V expansion for nonequilibrium quan- 
tum fields. We included contributions at N 2 LO, and 
found that the expansion is convergent even at large cou- 
pling as N is increased. At fixed N, we found convergence 
as the coupling strength is reduced. These results confirm 
that the 2PI-1/7V expansion truncated at NLO already 
gives quantitatively accurate results for moderate values 
of N or coupling strength: At any coupling N = 10 works 
very well, whereas for the often used N = 4, a trade-off 
in terms of a smaller coupling is required for precision 
calculations. We expect this behaviour to persist in 3 + 1 
dimensions, but a direct confirmation would be welcome. 

The numerical work was conducted on the Murska 
Cluster at the Finnish Center for Computational Sciences 
CSC. G.A. and N.L. are supported by STFC. A.T. is sup- 
ported by Academy of Finland Grant 114371. 
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